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Abstract 

We present elastic and inelastic spin-changing cross sections for cold and ultracold NH(X 3 E~) + 
NH(X 3 £~) collisions, obtained from full quantum scattering calculations on an accurate ab initio 
quintet potential-energy surface. Although we consider only collisions in zero field, we focus on the 
cross sections relevant for magnetic trapping experiments. It is shown that evaporative cooling of 
both fermionic 14 NH and bosonic 15 NH is likely to be successful for hyperfine states that allow for 
s-wave collisions. The calculated cross sections are very sensitive to the details of the interaction 
potential, due to the presence of (quasi-)bound state resonances. The remaining inaccuracy of 
the ab initio potential-energy surface therefore gives rise to an uncertainty in the numerical cross- 
section values. However, based on a sampling of the uncertainty range of the ab initio calculations, 
we conclude that the exact potential is likely to be such that the elastic-to-inelastic cross-section 
ratio is sufficiently large to achieve efficient evaporative cooling. This likelihood is only weakly 
dependent on the size of the channel basis set used in the scattering calculations. 
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I. INTRODUCTION 



Cold (T < 1 K) and ultracold (T < 1 mK) molecules offer a wide variety of applications 
in condensed-matter physics high-precision spectroscopy (2-4], physical chemistry [5 9 ] , 
and quantum computing [lOj, [ll|. In the last few years, techniques have been developed that 



either form (ultra) cold molecu 
112 1 or Feshbach association 



es by pairing up pre-cooled atoms, e.g. by photoassociation 
13] , or by cooling; the molecules directly. Examples of the latter 
approach include buffer-gas cooling [14| and Stark deceleration [15]. 

A promising candidate for direct-cooling experiments is the NH biradical. NH in its 
electronic X 3 X~ ground state has been cooled from room temperature using a helium buffer 



gas and trapped in a magnetic field jsl. [l6l 18j] . Stark deceleration and electrostatic trapping 
experiments have been performed on metastable NH(a 1 A), which, in contrast to the X 3 X~ 
ground state, exhibits a linear Stark effect. The decelerated NH(a 1 A) molecules can be 
converted to the ground state by excitation of the A 3 U a 1 A transition followed by 
spontaneous emission [l9|, |2(|. The resulting NH(X 3 £~) molecules may subsequently be 



accumulated in a magnetic trap. 

Direct-cooling techniques for NH are currently limited to temperatures of a few hundred 
mK. Producing NH molecules in the ultracold regime requires a second-stage cooling mech- 



anism, e.g. sympathetic cooling with an ultracold atomic gas |2lN25| or evaporative cooling. 
The latter process relies on elastic, thermalizing NH + NH collisions as the magnetic trap 
depth is gradually reduced. Inelastic spin-changing collisions lead to immediate trap loss 
and are therefore unfavorable. It is generally accepted that, in order to achieve evaporative 
cooling, elastic collisions should be much more efficient than inelastic ones. More specifi- 
cally, a Monte Carlo study on evaporative cooling of cesium atoms indicated that the ratio 
between elastic and inelastic collision rates should be greater than 150 26|. Although evap- 
orative cooling of NH might work with a lower ratio, it will be assumed that 150 is also the 
minimum required value for NH + NH collisions. 

For two magnetically trapped NH(X 3 E~) molecules, the collision complex is in the low- 
field seeking \S = 2,M S = 2) quintet spin state, with S denoting the total electronic spin 
and Ms its projection on the magnetic-field axis. Inelastic transitions may change either the 
Ms quantum number of the quintet state, or the total spin S to produce singlet or triplet 



complexes. The 5 = and 1 dimer spin states are chemically reactive 27H29I] and could be 
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of interest in cold controlled chemistry experiments (9]. 

A rigorous calculation of elastic and inelastic cross sections requires a full quantum 
coupled-channels method. In the case of NH-NH, however, the strong anisotropy of the 
interaction potentials and the open-shell nature of the monomers gives rise to a very large 
number of channels, making the calculation extremely challenging. In a recent study by 
Tscherbul et al. [30j on the iso-electronic 02(X 3 £~) - 02(X 3 £~) system, close-coupling 
calculations were performed that included up to 2526 channels, yielding cross sections con- 
verged to within 10%. These calculations were carried out in a fully decoupled channel basis 
to study collisions in the presence of an external magnetic field. It was noted, however, that 
the true O2-O2 interaction potential is likely to be more anisotropic than the potential used 
in their work, thus implying that even more channels would be needed. Other quantum 
scattering studies on O2-O2 include those by Avdeenkov and Bohn and Perez-Rfos et 



al. 



32| . In the work of Avdeenkov and Bohn, field-free collisions were studied using a total 



angular momentum representation, thereby reducing the total number of channels to 836. 
The rotational basis-set size used in these calculations was, however, smaller than that used 
in Ref. 



3CM ] . Perez- Rfos et al. also employed a total angular momentum basis, but the O2 
monomers were treated as closed-shell molecules. This allowed them to reduce the number 
of channels to 300. 

To our knowledge, only one theoretical study has been reported for the NH-NH system. 
Kajita SiJ employed the Born approximation, distorted-wave Born approximation, and 
classical path method to calculate elastic and inelastic cross sections at energies ranging 
from 1 /iK to 10 K, and found that evaporative cooling of NH is likely to be feasible. It 
must be noted, however, that only the electric dipole-dipole and magnetic dipole-dipole 
interactions were considered in these calculations. 

The aim of the present work is to obtain cold and ultracold NH + NH collision cross 
sections from rigorous quantum scattering calculations on an accurate ab initio quintet 
potential-energy surface. We include intramolecular spin-spin, spin-rotation, and intermolec- 
ular magnetic dipole-dipole coupling in the dynamics. In addition, we seek to address the 
issue of dealing with very large basis sets in order to converge the scattering results, a 
problem that is general for open-shell systems with relatively deep potential energy wells. 
For this purpose, we have employed a total angular momentum representation to perform 
the scattering calculations, assuming zero field. Collisions in a magnetic field are discussed 



3 



in a separate publication [34|. It will be shown that, within the uncertainty limits of the 
interaction potential, even an unconverged basis set can provide meaningful results. 

This paper is organized as follows. In Sec. [TTJ we discuss the scattering Hamiltonian and 
channel basis-set functions, followed by the details of the cross-section calculations. Results 
are presented in Sec. IIII Al In Sees. IIII Bl and IIII C\ we provide a comprehensive discussion 
on the accuracy of our calculated cross sections. Conclusive remarks are given in Sec. IIVI 

II. THEORY 

A. Hamiltonian and channel basis functions 

We consider the case of two colliding NH( 3 £~) molecules in the absence of an external field 
and treat the monomers as rigid rotors. We use a space-fixed coordinate frame to describe 
the collision complex. The relevant Jacobi coordinates are the intermolecular vector R that 
connects the centers of mass of molecules A and B, and the polar angles Ui = (6^,^) of 
the two monomers (i = A,B). We will neglect hyperfme coupling and assume that both 
monomers are in their nuclear-spin stretched states \I, Mi = I), with I = Jn + In denoting 
the maximum total nuclear spin and Mj its laboratory-frame projection. For fermionic 14 NH 
the maximum nuclear spin is I = 3/2 and for bosonic 15 NH we have 1 = 1. 

The scattering Hamiltonian for NH-NH can be written as 

h 2 d 2 L 2 

^ = ~^24iRdR?^ + 2 R 2 + UAj Ub ^ + ^ma g n.di P (-R, Sa, S b ) + H A + H B , (1) 

where \x is the reduced mass of the complex, R is the length of the vector R, L 2 is the 
angular momentum operator associated with rotation of R, Vs(R,uja,ujb) is the potential- 
energy surface for total spin S, V maugn ,di P (R, Sa, Sb) is the intermolecular magnetic dipole 
interaction between the two triplet spins, and Ha and Hb are the Hamiltonians of the 
individual monomers. The magnetic dipole term is given by 

2 

Wdip(«, S A , S b ) = -V6g 2 s ^ ^(-1) 9 C 2 ,_ 9 (Q)[^ ® S B ]f\ (2) 

i 

where gs ~ 2.0023 is the electron g-factor, /xb is the Bohr magneton, a is the fine-structure 
constant, G<i q is a Racah-normalized spherical harmonic, Q = (O, $) describes the orienta- 
tion of -R in the space-fixed frame, and the factor in square brackets represents the tensorial 



product of the monomer spin operators Sa and Sb- The monomer operators Hi each contain 
a rotation, spin-rotation, and intramolecular spin-spin term: 

Hi = B Q N? + 7iV, • + ~V6A SS ^(-l) 9 C 2 ,_ g (^)[.% ® Si]®, (3) 

g 

with Ni denoting the rotational angular momentum operator of monomer i. For brevity, 
we will denote the intramolecular spin-spin operator as Vgg . The numerical values for the 
rotational, spin- rotation, and spin-spin constants are B = 16.343275 cm -1 , 7 = —0.05486 
cm -1 , and Ass — 0.91989 cm -1 [35J for 14 NH, and, by scaling with the isotope mass (see e.g. 



p. 239 of Ref. [36|), we obtain B = 16.270340 cm" 1 , 7 = -0.05460 cm" 1 , and A ss = 0.91989 
cm" 1 for 15 NH. 

For the interaction potential Vs(R,oja,wb) we take the S = 2 ab initio surface of Ref. 



281 ] . This spin state corresponds to the case where both molecules are in their magnetically 
trapped (spin-stretched) states. Although the potential is based on the Jacobi coordinates 
for 14 NH - 14 NH, we use the same surface for the 15 NH - 15 NH isotope. This approximation 
is very reasonable since the center of mass of 15 NH is shifted by only 0.008 a$ with respect 
to that of 14 NH. We have verified that, at the equilibrium distance of the complex, this 
would give a maximum error of 2.2% in the 15 NH - 15 NH potential, which falls within the 



uncertainty range of the ab initio data. Following Ref. 



• 3, 



we expand the quintet potential 



in terms of spherical harmonics Yl,m °f degree L and order M: 

V(R,u a ,ujb)= Yl v l a ,l b ,l ab (R)A La , LbiLab (Q,uja,uJb), (4) 

La, Lb , Lab 

(L A MaL b Mb\LabMab) 

MaMbMab 

x YlaMa {ua)Yl b ,m b ( uj b)YI abMab (ft), (5) 

where (LaMaL b M b \LabMab) is a Clebsch-Gordan coefficient and the superscript * denotes 
complex conjugation. The subscript S = 2 has been omitted for brevity. It should be noted 



that the angular functions of Eq. (jSJ) differ by a factor o: 



C = ( _ 1) L A -L fl(47r) -3/2 (2LAs + 

1)[(2La+1)(2Lb + 1}] 1 / 2 from the functions used in Ref. 28[], i.e. the vla,l b ,l ab (R) expansion 
coefficients of Ref. [28] must be multiplied by ( to obtain the potential in the form of Eq. 

In the absence of an external field, both the total angular momentum J and its space- 
fixed projection Ai are rigorously conserved. We therefore expand the wave function in a 



total angular momentum basis: 

* J > M (R,n,UA,UB,<TA,<TB) = -B Y X J n^n b ,N,S a ,S b ,S,.J,l( R ) 



R 

N A ,N B ,N,S A ,S B ,S,J,L 



X 



' 1 1 J N a ,N b ,N,S a ,S b ,S,J,l(^', U A, UB, 0A, &B ), (6) 



where a a and Ob refer to the electronic spin coordinates of molecules A and B, respectively. 
Here Na and N B denote the rotational quantum numbers of the two monomers, N is the 
coupled rotational quantum number of the complex, Sa and Sb are the monomer spin 
quantum numbers, which are coupled into total spin S, J is the angular momentum quantum 
number arising from the coupling of N and S, and L denotes the partial-wave angular 
momentum. The coupled angular momentum basis functions are defined as follows: 

^N a ^ b ^,S a ,S b ,S,J,l(^ u a, u B , OA, ob) = 

Y Y Y Y y na,mn a ( uj a)y NbMnb (u b )y LMl (^) 

Mj,Ml M n ,Ms M Sa> Ms b M Na ,Mn b 

xt Sa ,m Sa {oa)t Sb ,m Sb (<?b) (N a M Na N b M Nb \NM n ) (S a M Sa SbM Sb \SM s ) 
x(NM n SM s \JMj)(JMjLM l \JM), (7) 

where ts a ,m Sa and ts Bi m Sb are spinor wave functions. Here the quantum numbers Mjv i; Ms t , 
Mat, Ms, Mj, and Ml denote the projections of Ni, Si, N, S, J, and L onto the magnetic- 
field axis. We will restrict the basis such that Na and iV^ range from to iV max and 
L = 0, . . . , L max . Note that the scattering calculations in this basis may also be performed 
for a single dimer spin state S. As detailed in Section IIIIBl we will exploit this feature 
to investigate the validity of describing all three dimer spin states by the 5 = 2 potential 
energy surface. 

Since target and projectile are identical, we can symmetrize the wave function with respect 
to the permutation operator Pab- This yields the following normalized basis functions: 

<Pn a ,n b ,n,s a ,s b ,s,j,l - [2(1 + SnaNbSsaSb)] i/2 [Vn a ,n b ,n,s a ,s b ,s,j,l 

+ ^(-l) i+ ^ +JVs - JV+SA+5fl - S ^;^,iV ! 5 fl ,5 A) 5, J) J • (8) 

Here rj — +1 corresponds to composite bosons and r] = —1 to composite fermions, assuming 
that the molecules are in their nuclear-spin stretched states. To obtain a linearly independent 



basis, the index pair (Na,N b ) must be restricted such that Na > N B |37|. Finally, the 



basis functions of Eq. (jSJ) are also eigenfunctions of the inversion operator, with eigenvalues 
e = (— \^ n a+Nb+l _ Thus, the Hamiltonian in the symmetry-adapted basis consists of four 
blocks, each block labeled by rj and the parity e. It must be noted, however, that the wave 
function of Eq. (JSJ) vanishes for (77 = +l,e = —1) and (77 = — 1, e = +1) if the molecules are 
in the magnetically trapped ground state with N A = N B = N = and S = 2. We therefore 
only need to consider the parity case e = +1 for 77 = +1 and e = — 1 for 77 = — 1. 

The matrix elements of the Hamiltonian in the symmetry-adapted basis [Eq. (jSJ)] can be 
readily obtained from the matrix elements in the 'primitive' basis iPn'^n b n s a s b s j l- These 
are given in the Appendix. 



B. S'-matrices and cross sections 



The close-coupling equations are solved for each J and each symmetry type (77, e) using 



the hybrid log-derivative method of Alexander and Manolopoulos [38] . This algorithm uses 
a fixed-step-size log-derivative propagator in the short range and a variable-step-size Airy 
propagator in the long range. The solutions are then matched to asymptotic boundary 
conditions to obtain the scattering S'-matrices. Since we consider only the field-free case, 
the results are independent of the total angular momentum projection Ai. 

Although we assume zero magnetic field in our calculations, we are ultimately interested 
in the elastic and inelastic spin-changing cross sections for magnetically trapped NH. It 
is therefore necessary to transform the S'-matrices to a channel product eigenbasis of the 
form \(Na, Sa)Ja, Mj a )\(Nb, Sb)Jb, Mj b )\L, Ml), where Jj and Mj i arise from the angular 
momentum coupling of Aj and SV Here we have used the notation A, instead of iVj, because 
Ni is strictly not a good quantum number. This is due to the intramolecular spin-spin 
coupling, which mixes states with Ni and Ni ±2. However, the mixing is quite weak and Ni 
corresponds almost exactly to Ni. A symmetry- adapted version of the channel eigenbasis is 
given by 

.„ 1 

^n a ,s a j a m Ja ,n b ,s b j b ,m Jb ,l,m l [2 (i + 6 NaNb 5 Sa sJj a j b 5 Mja m Jb )]^ 

x [\{N A , S a )Ja, M Ja )\{N b , S b )J b , M Jb )\L, M l ) 

x +r ] (-l) L \(NB,S B )JB,Mj B )\(N A ,S A )J A ,Mj A )\L,M L )}. (9) 
It should be noted that the total angular momentum J is not a good quantum number here, 



but its laboratory-frame projection M. = Mj A + Mj B + Ml is conserved. 

The basis transformation from Eq. (jSJ) to Eq. ([9]) cannot be performed analytically, be- 
cause Ni, N, and S are only approximately good quantum numbers. We have therefore 
developed a numerical scheme in which the channel eigenfunctions of Eq. are obtained 
as the simultaneous eigenvectors of the operators {L 2 , Ha + H B , J ZA + J ZB , J 2 A + J 2 }. Note 
that these operators all commute with each other and with Pab- The numerical procedure 
works as follows. We start by diagonalizing the first operator, e.g. the matrix representation 
of the L 2 operator, constructed in the basis of Eq. ((SJ). In each degenerate subspace of 
L 2 , we set up the matrix of the next operator and diagonalize it. This process is repeated 
for the remaining operators until all eigenvectors are unique. We note that the operator 
J 2 A + J 2 B is only required to distinguish between states with coincidental degeneracies in 

M Ja + M Jb , e.g. the states |^o,i,i,o,o,i,i,o,o,o> with m Ja = M J B = and I0o,i,i,i,o,i,i,-i,o,o) 
with Mj A = l,Mj B = —1. Any remaining degeneracies arising from Ha + Hg may be 
lifted by diagonalizing the operator H 2 A + H B , but such degeneracies occur only for higher 
energies. In the cold and ultracold regime, these higher-energy channels are closed and the 
eigenvalues of L 2 , Ha + Hb, J Za + Jz B , and J 2 + J 2 B are sufficient to uniquely identify 
all relevant quantum numbers. It must be noted that, since J ZA and J ZB do not separately 
commute with Pab, the matrices of J ZA + J ZB and j 2 + J 2 B are not trivially constructed in 
the basis of Eq. (JS}. We obtained these matrices by first evaluating the J Zi and J 2 , operators 
in a fully decoupled basis of the form \Na, M^ a , Sa, Ms a , Nb, Mn b , S b , Ms b , L, Ml). Both 
J Zi and J 2 , are diagonal in this basis, with diagonal elements Mj i = M^ + and Mj., 
respectively. We subsequently performed an analytical transformation to the coupled basis 
of Eq. flTJ) using the appropriate Clebsch-Gordan coefficients. Finally, we used a rectangular 
transformation matrix for J ZA + J ZB and J 2 A + J Zg to account for the symmetry adaptation, 
i.e. to transform the matrices to the basis of Eq. (jSJ). 

The evaporative cooling rate for cold magnetically trapped NH molecules, with quantum 
numbers Na = Nb = 0, J a = Jb = 1, and Mj A = Mj B = 1, is determined by the ratio 
between elastic and Mj-changing cross sections. The cross-section expression for indistin- 



guishable molecules at total energy E is 



30] 



a iAlB^l' A l' B ^ ~ ^ U2 7A7S ' ) \ T jAl B LM L ;y A j' B L'M' L ( E ) > ( 10 ) 

lAlB L ,M L L',M'r 



where we have introduced the shorthand notation 7,475 to label the symmetrized monomer 



states, i.e. (/>% A! s a ,j a ,m Ja ,n b ,s b ,j b ,Mj b ,l,m l = hAlB)\LM L ), and k lAlB is the length of the 
wavevector for the initial collision channel 17,47b)- The T-matrix elements are defined in 
terms of the transformed S'-matrix elements as T^ a ^ bLM ^, y UM , = ^-y A y A ^-y B -y' B ^ 'll> '<5 'm l m' l ~ 
^iai b lm l -i' i L'M' ■ F ma hy, we note that the summations over Ml and M' L in Eq. (HOP may 
also be understood as a sum over all possible M. values, since M. = M Ja + M Jb + M L = 
M' Ja + M' Jb +M' l . 



C. Computational details 

The scattering calculations were performed using a modified version of the MOLSCAT 



package [39j, |40| in which the coupled basis set of Eq. ([6]) was implemented. The radial 
grid ranged from 4.5 to 500 do, with the Airy propagation starting at 15 a^. The step size 
for the log-derivative propagator was 0.02 oq- The basis set included all functions up to 
Na = Nb = 5 and L = 6. The expansion of the quintet potential was truncated at La = 
Lb = 6. As mentioned in Section Hi Al the chemically reactive singlet and triplet interaction 
potentials were excluded from the calculations, and were replaced by the nonreactive S = 2 
surface. Thus, we assumed that all three spin states are described by the same potential 
energy surface. In order to study the role of the S = and 1 states under this assumption, 
we also performed scattering calculations for the quintet state only. 

At each collision energy, the scattering S'-matrices were accumulated for all relevant J 
values and subsequently transformed to the channel eigenbasis of Eq. d^J) for all possible 
M. values. The basis transformation was carried out in Matlab [41]. The total elastic and 
inelastic cross sections were then obtained using Eq. ([101) . 



III. RESULTS AND DISCUSSION 



A. Cross sections 

The elastic and Mj-changing cross sections for magnetically trapped 14 NH and 15 NH are 
shown in Fig. [TJ At low collision energies, the cross sections are dominated by incoming s- 
waves for bosonic 15 NH and by p-waves for fermionic 14 NH. The observed energy dependence 
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is consistent with Wigner's threshold law for iso-energetic processes 42J, |43) : 

aocE L+L ', 



(11) 



where L and L' denote the partial waves in the incoming and outgoing channels, respectively. 
For elastic 15 NH + 15 NH collisions, we have L = L' = and the cross section is constant as 
a function of E. For inelastic collisions, the change in Mj A or Mj B must be accompanied 
by a change in the Mi quantum number, which follows from the conservation of Ai. Since 
the parity (—1} n a+n b +l j g a j go r jg orous iy conserved, it is easily verified [see Eq. (jSJ)] that 
the dominant inelastic cross section for 15 NH {L = 0) corresponds to the V = 2 outgoing 
channel, and consequently behaves as E 2 . For fermionic 14 NH + 14 NH collisions, both the 
elastic and inelastic channels are dominated by L = L' = 1 [see Eq. (jSJ)], yielding the 
observed E 2 behaviour. We also point out that, in the presence of a magnetic field, all 
inelastic transitions would be exothermic and the corresponding cross section would behave 



as 



E L - 1 / 2 



42J. This leads to a different elastic-to-inelastic collision ratio than in the field- 
free case. It is shown in a separate publication that the ratio for 15 NH + 15 NH collisions is 



still very favorable when the magnetic field is explicitly included 



34] 



We find that NH is more suitable for evaporative cooling than NH, in agreement with 



the findings of Kajita 33]. More specifically, we see in Fig. [T] that the elastic-to-inelastic 
ratio for 15 NH + 15 NH far exceeds the critical value of 150 for all energies below E ~ 10 -2 
K, while for 14 NH + 14 NH the ratio is orders of magnitude smaller and is close to unity at 
collision energies below 10~ 4 K. This result is essentially a consequence of the Pauli principle, 
which forbids s-wave scattering for 14 NH + 14 NH. We emphasize that our calculations were 
performed under the assumption that both molecules are in their nuclear-spin stretched 
states, giving rise to a symmetric nuclear-spin wave function. This leads to the restriction 
that 7] = +1 (e = +1) for 15 NH and 77 = -1 (e = -1) for 14 NH. If, however, the two 
monomers were in different nuclear-spin states, the corresponding wave function may also 
be antisymmetric under exchange and both values of 77 would be allowed. In that case, the 
total cross section is given by a weighted sum over the cross sections cr +1 and a -1 : 

= W+ <L^* B W + W-<tb^m,W> ( 12 ) 

with W + and W~ denoting the relative spin-statistical weights. The weights are W + = 5/12 
and W~ = 7/12 for fermionic 14 NH and 3/4 and 1/4 for bosonic 15 NH. Figure |2] shows the 
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results for 14 NH - 14 NH, assuming a mixture of different nuclear-spin states. The inclusion 
of even-L partial waves (rj = +1) strongly enhances the efficiency of evaporative cooling for 
14 NH, in particular due to the s-wave elastic contribution. For 15 NH - 15 NH, the addition 
of odd-L partial wave contributions (rj = —1) will probably lead to a slightly lower elastic- 
to- inelastic ratio. This is because the odd-L elastic cross section, which vanishes as E 2 , is 
almost negligible compared to the s-wave elastic cross section in the ultracold limit. The 
odd-L inelastic contribution, on the other hand, exhibits the same threshold behaviour as 
the even-L inelastic cross section, and could easily increase the total inelastic loss by a factor 
of ~ 2. Hence we conclude that, in order to achieve efficient evaporative cooling, bosonic 
15 NH should be prepared in a single nuclear-spin state, while for 14 NH the molecules should 
be in a mixture of hyperfine states. 

Aside from symmetry arguments, the difference between 15 NH - 15 NH and 14 NH - 14 NH 
is relatively small. The rotational and spin-rotation constants differ by only 0.45% and the 
reduced masses of the collision complex are 6.6% different. Since 15 NH is more advanta- 
geous for evaporative cooling, we will only consider collisions between 15 NH molecules in 
the remainder of this work. Again it will be assumed that the monomers are in identical 
hyperfine states, so that only the rj — +1 (e = +l) symmetry case needs to be examined. 

State-to-state inelastic cross sections for magnetically trapped 15 NH (Mj A = 1, Mj B = 1) 
are shown in Fig. [3j We find that transitions to the states with \Mj A = l,Mj B = 0), 
\Mj A = l,Mj B = —1), and \Mj A = 0,Mj B = 0) are dominant in the ultracold regime. It 
can also be seen that these cross sections follow an E 2 dependence below ~ 10~ 4 K. The 
inelastic cross sections for \Mj A = 0,Mj B = —1) and \Mj A = —l,Mj B = —1) exhibit E 4 
behaviour at low collision energies. These results are consistent with the threshold laws of 
Krems and Dalgarno for collisional reorientation of angular momentum in the absence of an 



external field [43j. Although these laws were derived for collisions of paramagnetic species 
with structureless targets, they also apply to 15 NH + 15 NH collisions: 

<J J,Mj^J,Mj±AM J OC E AMj (13) 

if AMj is even and 

aj, Mj ^j,Mj±AM J cc E AMj+1 (14) 

if AM j is odd. Here AMj is defined as the change in Mj A + Mj B . It also follows from Eq. 
(fl3"|) that the elastic cross section (AM; = 0) is constant at low energies, in agreement with 
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Eq. flm>. 



B. Contributions from singlet and triplet states 

Throughout this paper, we have assumed that all three spin states of the NH-NH complex 
are described by a single nonreactive potential-energy surface, namely the S = 2 surface. 
The S — 2 state corresponds to the case where both monomers are magnetically trapped, 
and is therefore the most relevant spin state in our present study. It is, however, not a priori 
clear how the S = and 1 states can influence the trap loss probability, and how well they 
can be described by the quintet surface. 

We must first point out that, even at infinite separation, S is strictly not a good quan- 
tum number due to the intramolecular spin-spin coupling. However, the coupling between 
different spin states is relatively weak and we may therefore treat S as nearly exact. Specif- 
ically, for the rotational ground state of the complex, the initial state with Mj A = Mj B = 1 
corresponds almost exclusively (99.98%) to the quintet state. 

In order to investigate the contributions from the S = and 1 states, we have performed 
scattering calculations with all singlet and triplet functions removed from the basis set. The 
results are shown in Fig. H] as a function of energy. The cross sections for the full basis set, 
i.e. with all three spin states included, are also plotted for comparison. It can be seen that 
exclusion of the S = and 1 states has a rather small effect on the cross section, suggesting 
that most of the trap loss takes place within the quintet state. Thus, the singlet and triplet 
states play a minor role in the collision dynamics when described by the nonreactive S — 2 
potential. 

If the S = and 1 states would be described by their true, reactive surfaces, it can be 
expected that any transition to the singlet or triplet state leads to chemical reaction and 
consequent trap loss. In that case, however, the potentials are no longer degenerate at short 
range and the probability for hopping from the quintet surface to another state is most 
probably decreased due to the energy gap law. That is, inclusion of the reactive S = and 
1 surfaces will probably not lead to a larger inelastic cross section, and our assumption of 
including only the nonreactive S = 2 surface is very reasonable. In this respect, we may 
also view the Mj-changing cross sections presented in Fig. H] as approximate upper bounds. 
Nevertheless, it must be noted that the relatively deep wells in the reactive potentials will 
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give rise to a large number of bound states, which in turn may cause strong resonances in the 
cross sections. In order to verify these assumptions, we plan to perform reactive quantum 
scattering calculations for NH + NH with all three interaction potentials included. 



C. Sensitivity to potential and basis-set size 

In this section we address two interrelated topics, namely the sensitivity to the potential 
and the dependence on the angular basis-set size. It is well established that low-energy 
scattering depends strongly on the presence of bound and quasi-bound states near the dis- 
sociation threshold. Such states can give rise to scattering resonances that may enhance the 
collision cross section by several orders of magnitude. The energies of these (quasi-)bound 
states are highly sensitive to the details of the potential-energy surface, and hence they 
are very difficult to predict from first principles. Even a state-of-the art ab initio potential 
cannot reliably predict whether a particular near-dissociation state lies above or below the 
threshold. This is particularly true for systems with multiple degrees of freedom and deep 
potential wells, for which the density of states is relatively high. Thus, in order to assess the 
accuracy of the cross sections, we must carefully take into account the effect of uncertainties 
in the potential. In a related manner, we also consider the effect of using different channel 
basis-set sizes in the scattering calculations. The size of the angular basis set can influence 
the energies of the (quasi-)bound states, which in turn can lead to a different resonance 
structure. It will be demonstrated, however, that the use of a reduced basis set leads only 
to a shift in the resonance positions, and does not significantly alter the general resonance 
pattern. 

We first consider the sensitivity of the calculated cross sections to the potential-energy 
surface. Our potential has been obtained from state-of-the-art ab initio calculations, and 
we estimate that it differs from the exact potential by at most a few percent. For practical 
reasons, we have studied the potential dependence indirectly by performing scattering cal- 
culations as a function of the reduced mass /i. Since scaling the reduced mass by a factor 



of A (/i S caied = A/i) is almost equivalent to scaling the entire interaction potential by A |24(, 
this provides a stringent test for the sensitivity to the potential. The true potential does 
not necessarily differ from our ab initio surface by only a constant factor, but scaling by A 
(0.9 < A < 1.1) amply samples the range of possibilities within which the exact potential is 
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expected to lie. 

Figure [5] shows the cross sections as a function of A at collision energies of 10 -6 K, 10 -4 K, 
and 10 -3 K. It can be seen that both the elastic and inelastic cross sections change by several 
orders of magnitude as a function of A, but they vary about a certain background value. 
For instance, the elastic cross sections fluctuate around ~ 10~ 12 cm 2 for all three collision 
energies. The background values for the inelastic cross sections increase with E 2 in the 
ultracold regime, consistent with the results of Fig. [TJand the threshold laws discussed in Sec. 
IIII A[ The deviations from the background values are due to scattering resonances, which 
arise from NH-NH states that change from bound to quasi-bound at the \Mj A = 1, Mj B = 1) 
threshold. Such resonance features are to be expected as a function of A, since a scaling of 
the potential, or in fact any modification of the potential-energy surface, will cause a shift in 
the bound-state energies. For 10~ 6 , 10~ 4 K, and 10~ 3 K, the resonances are located around 
the same values of A, and hence the A-dependent resonance structure would not be averaged 
out in a thermal (Maxwell-Boltzmann) distribution at temperatures below 1 mK. That is, 
thermally averaged rate constants are likely to show a similar sensitivity to the potential as 
the calculated cross sections. 

Let us now consider the elastic-to-inelastic cross-section ratios as a function of A. These 
are shown in Fig. [6] for E = 10~ 6 K, 10~ 4 K, and 10~ 3 K. For clarity, we have also indicated 
the critical ratio of 150 that is required for efficient evaporative cooling. As can be seen, the 
calculated ratios exceed 150 for almost all values of A and all energies considered, except 
when A is close to resonance. This demonstrates that evaporative cooling of NH is feasible at 
energies below 1 mK for most of the A- values considered. Although we cannot predict which 
value of A corresponds most closely to the exact potential, we do expect that the sampled 
range of A is indicative of the range within which the exact potential lies, and hence we 
conclude that the probability for successful evaporative cooling is relatively large. That is, 
the true potential is very likely to be such that the elastic-to-inelastic ratio exceeds 150. 

The A-scaling approach is also used to investigate the influence of the angular basis-set size 
on the scattering results. First we point out that the strong anisotropy of the potential and 
the large reduced mass of NH-NH require relatively high values of the basis-set parameters 
Af max and L max . In addition, the triplet spins on the monomers increase the channel basis- 
set size by a factor of 9, making it highly challenging to achieve full basis-set convergence. 
Figure [7] shows the cross sections as a function of A for different values of A" max and L max at 
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a collision energy of 10~ 6 K. The maximum number of channels in these calculations ranged 
from 937 for iV max = 4 and L max = 6 {J = 4) up to 2382 for iV max = 6 and L max = 6 
{JJ = 5). It can be seen that the cross sections all vary by several orders of magnitude as a 
function of A, and for a given value of A the four basis sets can yield very different numerical 
results. However, the different cross sections vary about the same background values and 
the resonant features have similar widths for all four basis sets. Thus, a change in iV max or 
L max may cause a shift in the positions of the resonances, but the overall pattern is virtually 
unaffected. The estimated probability for successful evaporative cooling, i.e. the probability 
that the exact potential is such that the elastic-to-inelastic cross-section ratio exceeds 150, 
is therefore similar for all four basis sets. This can also be understood by considering that 
a change in the basis set only shifts the bound-state energy levels, similar to the effect of 
scaling the potential. 

The results of Fig. [7] demonstrate that the cross sections are almost, but not fully con- 
verged with respect to iV max and L max . Using a larger basis set is infeasible at present given 
the available computer power. A larger basis set would also require additional terms in the 
expansion of the potential anisotropy [Eq. (J4j)], making the calculation prohibitively expen- 
sive. Moreover, taking into account the uncertainty in the potential, even a fully converged 
basis set would not give really reliable numerical values due to the presence of (quasi-)bound 
state resonances. Since the exact form of the potential, and thus the precise locations of 
the resonances, are still unknown, the calculated cross sections are subject to an inherent 
degree of uncertainty that cannot be reduced by the use of a fully converged basis set. In 
this sense, full basis-set convergence will not necessarily yield a more accurate prediction of 
the true cross sections. On the other hand, the probability for successful evaporative cooling 
can be reliably predicted using an incompletely converged basis set, and hence we conclude 
that, even if full basis-set convergence could be achieved, this would not significantly alter 
our main qualitative results. We emphasize, however, that it is crucial to test the sensitivity 
to the potential in order to assess the accuracy of the calculated cross sections. As a final 
point, we note that the uncertainty limits of the potential could, in principle, be greatly 
reduced by measuring the cross sections experimentally. 
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IV. CONCLUSIONS 



We have carried out elastic and inelastic quantum scattering calculations on a state- 
of-the-art ab initio potential to study field-free NH + NH collisions at low and ultralow 
temperatures. The results indicate that, when the molecules are prepared in their nuclear 
spin-stretched states, bosonic 15 NH is more suitable for evaporative cooling than fermionic 
14 NH. This is a direct consequence of the Pauli principle, which forbids s-wave scattering for 
two identical fermions. The 14 NH isotope may also be successfully cooled, however, when 
the monomers are in a mixture of different nuclear spin states. 

We have assumed that all three spin states of the NH-NH complex are described by 
the nonreactive quintet surface. This approximation is shown to be reasonable, although a 
full reactive scattering calculation would be required to investigate the precise role of the 
chemically active singlet and triplet states. 

The collision cross sections are sensitive to the details of the interaction potential, because 
of the presence of quasi-bound states that cause scattering resonances. Since the exact 
interaction potential is unknown, this gives rise to a degree of uncertainty in the numerical 
cross sections. However, a sampling of the range of possibilities indicates that the exact 
potential is very likely to be such that the elastic-to-inelastic cross-section ratio is favorable 
for evaporative cooling. This result is only weakly dependent on the size of the channel basis 
set. In particular, the effect of using a reduced basis set is very similar to a scaling of the 
potential within its uncertainty. We conclude that even without full basis-set convergence, 
which is extremely difficult to achieve for systems such as NH-NH, we can provide valuable 
insight into the feasibility of evaporative cooling. This also offers hope for the theoretical 
treatment of other challenging open-shell molecule + molecule systems. 
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Appendix: Matrix elements 



In this Appendix, we present the matrix elements of the scattering Hamiltonian in the 
'primitive' basis ' 1 Pn a n b ns a s b sjl- The matrix elements in the symmetry-adapted basis 
can be obtained using Eq. For the angular functions of the potential we have 



i,J,M 



J,M 

N A ,N B) N,SA,S B ,S ) J,L\ fJi -L A ,L B ,L AB \r N ^ N i B ^ N ^ SA ^ B ^^j^ u 



&<?.<?' I — 

4:71 



3/2 



I At a . r, a . r, a cj I VW' 

A 

_^\N A +N B +N+S+L AB +J 



3'J'L') 



x [Lab] \ [L A ] [L b ] [N a ] [N' a ] [N b ] [N' b ] [N] [N 1 ] [L] [U] [J] [J'] 



x 



X 




N A N' a L a 
N b N' b L b 



N N' L 



AB 



(Al) 



with the factors in large round brackets denoting Wigner 3j symbols, the factors in curly 
brackets denoting 6j and 9j symbols, and [Q] = (2Q + 1). The intermolecular magnetic 
dipole term is given by: 



IK, 



J,M ly \ i J ,Jvi 

N A ,N B ,N,S A ,S B ,S,J,L\ v ma.ga.dip\YN' A ,N^,N',S A ,S B ,S' ! J',L' 

a 2 



3'J'L't 



- 5n a n a 5n b n> b 6 NN >VWghl^(-l) N+s ' +J+J ' +J 

x y/S A (S A + 1)S B (S B + 1)[S a ][Sb][S][S%J][J'][L][L'] 

Sa Sa 1 

xi I < > < > < S B S B 1 

S S' 2 




(A.2) 



The rotation operators for the two monomers (i = A, B) are completely diagonal in the 
angular basis: 



$NAN' A 6N B N' B 8NN'$SS'6jJ>5LL'BoNi(Ni + 1). 



(A3) 
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For the spin-rotation coupling terms we find 

(^n'^Nb^^Sb&J^INa " SA^W^N' B ,N',S A ,S B ,S',J',L') = 
A A A A ru( 1 \N A +N B +S A +S B +S+S'+J 

x s/N A {N A + 1)S A (S A + 1)[N A ][S A ][N][N'}[S}[S'} 

(n a n a i Us A s A 1 jj^i 

[ TV TV' 7V B J [ 5 5' 5 B J [ 5' 5 J J 

{^N a ,N b ,N,S a ,S b ,S,J,l\1 N B ■ SB\i>N ] ^N' B ,N>,S A ,S B ,S',J>,L>) = 

x a/ Nb(Nb + 1)5b(Sb + 1)[^b][S S ][AW[S][5'] 

x { N> N 1 1 { S ' S 1 1 { N N> 1 (A 5) 

X [ iV B jV B iV A J I 5b 5b 5 a J \ S' S J 

and, finally, for the intramolecular spin-spin operators Vgg we have 

Wn a ,N b ,N,S a ,S b ,S,J,L\ v SS Wn' a ,N' b ,N',S a ,S b ,S',J',L'I - 
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E(K) 

FIG. 1: Elastic and inelastic Mj-changing cross sections for 14 NH + 14 NH and 15 NH + 15 NH 
collisions, assuming that all molecules are in their magnetically trappable and nuclear-spin stretched 
state. 
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FIG. 4: Elastic and inelastic Mj-changing cross sections for magnetically trapped 15 NH obtained 
from scattering calculations with only the quintet state included in the basis. The cross sections 
calculated with all three spin states included ( "all S") are shown for comparison. 
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FIG. 5: Elastic and inelastic Mj-changing cross sections for magnetically trapped NH as a 
function of the scaling parameter A, calculated at collision energies of 10 -6 K, 10 -4 K, and 10 -3 
K. The elastic cross sections for 10~ 4 K are the same as for 10 -6 K. 
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FIG. 6: Elastic-to-inelastic cross-section ratios for magnetically trapped NH as a function of the 
scaling parameter A, calculated at collision energies of 10" 6 K, 10 -4 K, and 10 -3 K. The horizontal 
black line indicates the critical value of 150 that is required for efficient evaporative cooling. 
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FIG. 7: Elastic and inelastic Mj-changing cross sections for magnetically trapped 15 NH as a 
function of A, calculated for different basis sets at a collision energy of 10 -6 K. Solid lines correspond 
to elastic cross sections and dashed lines to inelastic cross sections. Different colors represent 
different basis sets. 
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